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Abstract 

Linear optical networks are devices that turn classical incident modes by a 
linear transformation into outgoing ones. In general, the quantum version of 
such transformations may mix annihilation and creation operators. We derive 
a simple formula for the effective Hamiltonian of a general linear quantum 
network, if such a Hamiltonian exists. Otherwise we show how the scattering 
matrix of the network is decomposed into a product of three matrices that 
can be generated by Hamiltonians. 
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Simple optical instruments Ij such as beam splitters or parametric amplifiers 
are characterized by linear input-output relations. The beam splitter transforms 
the annihilation operator of the incident light modes according to the classical laws 
of optical interference, i.e., by a linear transformation. The parametric amplifier acts 
like a phase-conjugating mirror, combining the annihilation operator of one incident 
mode with the creation operator of the other. Complex optical networks can be 
constructed from beam-splitters, mirrors and active elements such as parametric 
amplifiers 0EH3- 

Networks are essential to the optical communication of information and for novel 
ways of quantum- information processing E] . Furthermore, linear networks may 
possess interesting quantum-statistical properties when they are large, for instance 
as representations of the Ising model Ej and as examples of quantum localization 
[Oj. Here we derive a simple formula for the Hamiltonian of an arbitrarily large 
quantum-optical network, if such a Hamiltonian exists. Otherwise we show how the 
scattering matrix of the network is decomposed into a product of three matrices that 
can be generated by Hamiltonians. Our results allows to predict how the quantum 
state of the incident light modes is processed. Our theory contains as special cases 
the previously studied Hamiltonians of symmetric passive networks j3j, the theory 
of the beam splitter PUJ EH E2| and of the parametric amplifier [T?^ I14j. 

Consider n incident modes with the annihilation operators a& and the creation 
operators aj., the index k running from 1 to n. Suppose that the optical network 
produces n' outgoing modes with the corresponding operators a' k and a£. (As we 
shall see in (J7J) below, n' must equal n.) 

Canonical quantization of the linear input-output behavior of a classical linear 
network leads to the requirement that the mode operators are related to each other 
by the linear transformation 



(In this formula, a,k is representative for the vector formed by all annihilation op- 
erators, and similar for the other operators indexed by k.) S denotes the classical 
scattering matrix of the network; it has 2n' rows and 2n columns and must be of 
the form 



where A denotes the (untransposed) complex conjugate of A. Indeed, (j2J) is necessary 
for the consistency of (Q), as can be seen by writing in components and forming the 
adjoint. If B = 0, the network is passive and can be constructed from beam splitters, 
phase shifters and mirrors [2]. Otherwise, the transformation ([IJ mixes annihilation 
and creation operators, and we speak of an active network. 

In the Schrodinger picture, incident light in the state l^) is transformed into the 
state l^') = S\i[>) of the outgoing light, with a unitary scattering operator S. In the 
Heisenberg picture, this corresponds to the transformation O' = S^OS for arbitrary 
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operators O. In particular, to match the mode transformation we need to have 
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S. 



(3) 



Here it is understood that S acts on each creation and annihilation operator sep- 
arately. In the following we show how to construct from the classical scattering 
matrix S a unitary scattering operator S such that Q holds; this operator then 
completely specifies the desired quantum behavior of the network. 

The light of both the incident and the outgoing modes consists of bosons subject 
to the commutation rules 
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We write the commutation relations (HI) in matrix form, 
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where 1 and 1' denote the unit matrix in n and n' dimensions, respectively. Inserting 
the linear mode transformation (0) into the commutation relations and using 
we find that the commutation relations are equivalent to G' = SGSt This relation 
implies that GS^G' is the inverse of S. Since S is invertible, the number of incident 
modes must equal the number of outgoing modes, 



n = n 



(7) 



a result well-known for passive networks made of beam splitters and mirrors 01 
Ej. This feature therefore remains true for general quantum-optical networks. In 
particular, even if just one light beam is transformed into two or more beams a 
matching number of incident vacuum modes are involved with their vacuum noise 
affecting the outgoing light. 



Since (jJJ) implies G' = G, we find 



G = SGS 1 



as another consistency condition for classical scattering matrices. The matrices S 
satisfying (jHj) are called quasi-unitary ^3] and form a Lie group Go whose infinitesi- 
mal generators are the elements of the quasi-unitary Lie algebra L consisting of all 
matrices K with 

KG + GK f = 0. (9) 
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Because of (JHJ), (jUJ) implies that the elements of L are precisely those of the form 



with antihermitian A, B and arbitrary D. Because of (J2J), an S suitable as a classical 
scattering matrix in fact belongs to a subgroup G of Go; the corresponding Lie 
subalgebra L consists of all K of the form 

with antihermitian A and complex symmetric D. 

By standard results for Lie groups, every element S G G can be written (in many 
ways) as a product 

S = e Kl • ■ • e K '" (12) 

of finitely many exponentials of infinitesimal generators K 1; . . . , K m G L. 

A special case in which a single exponential often suffices is when S G G can be 
diagonalized, i.e., S = XAX^ 1 with a diagonal matrix A, and has no eigenvalues 
that are real and negative. In this case S = e K with K = X(\nA)X~ 1 , where, for 
the logarithm In A, the principal value is taken in each diagonal element. Often, and 
if K is sufficiently small, always, K G L; then Q holds. 

For general S G G, a decomposition (fT2"|) of the form 

At \ / D\ fA 3 
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with antihermitian Aj and real diagonal D can be found constructively (see Ap- 
pendix A); the resulting Kj belong to L since they have the form (JUJ. The fac- 
torization can be given a natural physical interpretation, similar to the one of the 
known factorization of two- mode parametric amplifiers PjJ El El : Any linear net- 
work can be thought of acting in three steps. First the incident modes are mixed by 
a passive network, then they undergo parametric amplification with real squeezing 
parameters, and finally they are subject to another passive mode transformation. 
In the special case when the entire network is passive the squeezing parameters are 
zero. 

Given the structure of general group elements S, we may concentrate in the 
construction of the unitary scattering operator S satisfying © on the case of a 
single exponential S = e K with K G L. We observe that the operator 

£ K = ^(4A) H (^t ) withH = -*GK (14) 

is Hermitian. Indeed, H is Hermitian since 

H f = iK+G 1 = iK^G = -iGK = H (15) 
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by and the hermiticity of follows. The operator plays the role of an 
effective Hamiltonian in a fictitious time evolution for the quantum-optical network 
generating the linear mode transformation in the Heisenberg picture. (In 
simple cases, the fictitious time r corresponds to the fractional depth of the optical 
device.) To show this we put 



S(t) = exp(— itHk) 
and verify that the operator vectors 
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satisfy the same differential equation 

d 



dr 



A(t) = KA(t). 



(16) 
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Indeed, (fTSjl holds trivially for A(r). To show that it also holds for A{r), we use the 
commutation relations in matrix form © and get, using (fT7)|) . 
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so that the differential equation (|18|) also holds with A(t) in place of A{t). Since 
A(t) and A{j) trivially agree at r = 0, they agree for all r. In particular, we 
conclude that A(l) = A(l). In view of ()16|) and ()17j) . this implies that the unitary 
operator 

5 K = exp(-i^ K ) (20) 

satisfies the required equation (jSJ) for S = e K . Therefore, Sk is the desired scattering 
operator corresponding to the classical scattering matrix S = e K . 

In the more general case where the classical scattering matrix S is given by a 
product (fT2j) of exponentials, one sees by direct substitution that (jHJ) is satisfied by 
the scattering operator 

S = S Km ---S Kl . (21) 

Our general formula contains as special cases the known effective Hamiltonians 
for beam splitters and parametric amplifiers [1, Section 3.3] that are characterized 
by the scattering matrices 

/ cos (j) — sin (ft \ 
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with real (p for a beam splitter and 
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(23) 



with real ( for a parametric amplifier. One can easily diagonalize the matrices 
S S put (<j>) and S amp (() and calculate so their logarithms K. sp ut and K amp . However, 
the product 

S = Samp(OSspiit((fi) for cosh C cos = 1 (24) 

does not possess a diagonal representation. This case represents a beam splitter 
that is exactly compensated by parametric amplification. For example it describes 
an eavesdropping attempt where quantum information is tapped by beam splitting 
followed by amplification to restore the light intensity. To find the Hamiltonian for 
this device, one has to resort to a Jordan decomposition of the scattering matrix. 
The resulting matrix logarithm is K = S — 1. Finally, we mention a simple example 
where the logarithm of the scattering matrix fails to satisfy the Lie condition Q 
that is essential for the hermiticity of the Hamiltonian (fT4*j) . The example is a single- 
mode squeezer followed by a 7r phase shifter, characterized by the scattering matrix 



/ cosh C sinh A . , Q 
ysinh C cosh (y ' 



(25) 



Since the trace of S is —2 cosh £ < —2, the matrix K = log S must fail to satisfy 
(JHJ) for arbitrary choices of the branch of the logarithms. Indeed, if © holds, K 
has the form (fTT]) with purely imaginary A, hence K has two real or two purely 
imaginary eigenvalues A1.2. Thus the trace of S = e K is e Al + e A2 , which cannot be 
a real number < —2. Consequently, there is no single Hamiltonian that generates a 
single-mode squeezer followed by a it phase shifter, although both subdevices possess 
effective Hamiltonians. 

To summarize, we have developed an explicit procedure how to calculate the 
Hamiltonian of a quantum-optical linear network. If the matrix logarithm K of 
the scattering matrix S satisfies the Lie condition (jHJ), the Hamiltonian is given by 
(JHJ). One can calculate K by diagonalizing S or, if this is not possible, using the 
Jordan decomposition of S. We have shown that the scattering matrix S of any 
linear quantum-optical network can be decomposed into three factors that can be 
diagonalized or are already diagonal. Note that the Hamiltonian (|14|) is not unique, 
because the matrix logarithm is multivalued. Furthermore, since any decomposition 
(|T2~j) leads to a realization of the network, there are many equivalent ways to design 
an optical network with a particular input-output relation. There are also many 
ways to assemble it from the basic building blocks [2J El H] > from beam splitters and 
parametric amplifiers. 
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Appendix A 

In this appendix we prove that every S G Go can be written in the form (J 12)) with 
three factors of the form 

K *=(J o). K '=(t < 26 > 

with antihermitian Ajf. and real nonnegative diagonal D. In case that S G G, the 
construction can be modified such that ffTTTf) hold with antihermitian Aj and real 
diagonal D. 

The proof is constructive and begins by partitioning the matrix S into four n x n 
submatrices, 

'o 1 o 12 V (27) 
The condition (JBJ) which expresses the fact that S G Go implies the equations 

SiiSii — S12S12 — 1) (2£ 



•Sii^i — S12S22, (29) 

^i^i — S22S22 = —1- (30) 

(J2EJ) implies that WS^xW 2 = x^SnS^x = x^S u S\ 2 x + x ] x = \\S\ 2 x\\ 2 + \\x\\ 2 > 
if x 7^ 0, and therefore that Su is invertible. It is always possible to factor S12 
into a product Su = UiSV^ (singular value decomposition; see, e.g., P2>|) consisting 
of unitary matrices U\, V2 and a nonnegative real diagonal matrix S. The matrix 
C := (S 2 + l) 1 / 2 is a real, nonnegative invertible diagonal matrix commuting with 
S and satisfying 

C 2 - S 2 = 1. (31) 

We now form the matrices 

Vx := S^U X C, U 2 := S22V2C- 1 . (32) 

This immediately gives Su = U\CVi and 622 = U2CV2 , moreover (|2^|) implies that 
S21 = U2SV1 . Inserting this into (}28|) and (J30j) shows that V\ and U2 must be 
unitary, and insertion into (|27j) gives 

„ _ (u x cvt \J X SV\\ _(u, o \ (c s\ M \ 

b " [uisvi u 2 cv}) " V u 2 J\s c){o vi) ■ m 
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Using (jHJ), it is easily seen that each factor in this factorization belongs to Gq. Due 
to their special form, they can be easily be brought into exponential form. Indeed, 
since unitary matrices U are normal and have eigenvalues of absolute value one 
only, they have a spectral factorization U = Q exp(i$>)Q> with unitary Q and real 
diagonal $, so that U = e with antihermitian K = iQQQ^ . Therefore we can 
find antihermitian Ajk such that Uk = e Alk and V k = e Ask for k — 1,3. If we also 
define the real diagonal matrix D := log(C + S) with componentwise logarithms 
(nonnegative since C > 1, S > 0) on the diagonal, it is easily seen that 

s =™v( A o X 2 ) exp (^ o) exp (1r a° 32 )' ^ 

with antihermitian Ajk, as asserted. 

Now suppose that S G G. Comparing (J2J) and (|3^j) . we find that we must have 

U^vl = U 2 CV}, UiSV\ = U 2 SVl (35) 

If the diagonal entries of C and S are all distinct, the singular value decomposition is 
known to be unique up to a diagonal matrix of phases. Therefore, there are diagonal 
matrices Qj with Q^Qj = 1 such that 

U 1 = U 2 Q 1 , QxV{ = Vj, Ui = U 2 Q 2 , Q 2 V\ = V}. (36) 

Clearly, this implies that Q := Q\ = Q 2 is real diagonal with Q 2 = 1; in particular, 
Q has diagonal entries ±1. Therefore 

/ UxCV} U X SQV[ \_fUiO\fC QS\M \ 
\7JiQSVt U lQ CQVl) V0 UJ{QS c)[o v\) ' W 

As before, we can find antihermitian Aj such that U\ = e Al and = e As . If we also 
define the real diagonal matrix D := log(C + QS) with componentwise logarithms 
on the diagonal, it is easily seen that (|13|) holds with antihermitian Aj and real 
diagonal D. 

Finally, if some diagonal entries of C or S coincide, we can perturb S slightly to 
remove the degeneracy, and obtain the same decomposition by a limiting argument. 

Note added: An equivalent factorization came to our attention after finishing 
the paper. 
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